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to occur through a spinodal instability from a super-cooled initial state. We consider the 
creation of net Chern-Simons number, which through the axial anomaly is equivalent to a 
^«\j . baryon asymmetry. CP violation is explicit in the scalar potential, but only in combination 

^.^ , with P-violation is it possible for an asymmetry to be generated. This is introduced through 
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^r^ I out the fermions in the theory. We perform real-time lattice simulations of the transition, 

^^ ' and find the coefficient of this term required for successful baryogenesis. 
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1 Introduction 

The origin of the baryon asymmetry in the Universe is a major unresolved issue in cos- 
mology, and constitutes together with dark matter, cosmological inflation and dark energy 
the only direct evidence of physics beyond the Standard Model. Baryon number violating 
processes do exist within the electroweak sector of the Standard Model, and this realiza- 
tion has prompted extensive development of the Electroweak Baryogenesis scenario [1]. 
Within the Standard Model, one encounters two stumbling blocks to successful baryogen- 
esis: CP-violation as encoded in the CKM quark mixing matrix is much too small at high 
temperatures (T ~ m^) [2] and the electroweak phase transition is too smooth to supply 
the out-of-equilibrium conditions necessary for an asymmetry to be created [3]. 

Three main avenues are pursued to remedy these shortcomings. Firstly, a non-trivial 
mass sector for the leptons may provide for leptogenesis [4, 5], with additional CP-violation 
in the lepton mass matrix and heavy masses allowing for out-of-equilibrium decay. Sec- 
ondly, one may extend the Higgs sector to allow for a strongly first-order electroweak 
transition as well as CP-violation through complex couplings, reviving the original "hot" 
electroweak baryogenesis mechanism, see for instance [6, 7]. And thirdly, additional scalars 
may allow for the Universe to super-cool below the electroweak scale, and then trigger 
a low-temperature spinodal transition. Over the last decade, this has become known as 
"cold" electroweak baryogenesis [8-11]. 

The "cold" scenario was originally an attempt to bypass the lack of a strong equilib- 
rium phase transition by invoking an out-of-equilibrium process such as preheating after 
inflation, and necessary CP-violation was introduced through a hypothetical higher order 
bosonic operator [8, 11]. Through extensive numerical simulations, it was shown that this 
could lead to successful baryogenesis [11-13]. 



Later, Standard Model CP-violation itself was reintroduced in this "cold" setting. As 
was argued in [14-17], and more explicitly computed recently [18-20], the effect is not 
strongly suppressed at low temperatures. It was argued that if the transition has an 
effective temperature around 1 GeV, Standard Model CP-violation may be able to provide 
for baryogenesis. This awaits further study. 

In the 2HDM, which is our main concern in the present work, we will neglect the 
CKM matrix. Instead, CP-violation is present in the scalar interactions. In addition, 
the potential has to be such that the potential exhibits a strong spinodal instability (see 
below). We will assume that the Higgs fields are coupled to the Standard Model fermions, 
but include these only through what we expect to be the leading (C- and P-breaking) 
bosonic term arising from integrating them out of the path integral. The reason for this 
is that for bosonic simulations, the asymmetry is generated as non-zero average Chern- 
Simons number, which is both CP- and P-odd. The scalar potential breaks CP through 
breaking C, and so P is conserved. So we need the effect of the explicit P-breaking present 
in the electroweak gauge- fermion interactions, in our case through a bosonic C- and P- 
breaking operator^. 

Our goals are therefore threefold: To simulate a cold tachyonic symmetry breaking 
transition in the 2IIDM. To study the dynamical interplay between a source of C-violation 
and one of C- and P-violation to generate effective P- and CP-violation, necessary for the 
creation of a baryon asymmetry. And finally to explicitly compute the magnitude of the 
asymmetry as a function of in particular the coefficient of the C- and P-breaking term. 

Electroweak baryon number creation: In the electroweak sector of the Standard 
Model, an anomaly relates the fermion (baryon B, lepton L) numbers to the evolution of 
the Chern-Simons number Ncs of the SU(2) gauge field 

B{t) - B{0) = L{t) - L(0) = 3 [iVes(t) - iVcs(O)] =^ j dt j d'x ^TrF^,F^^ (1.1) 

with F^^ = \l2e^^^" Fpfj. The question of baryogenesis then becomes a question of whether 
permanent change of Chern-Simons number can take place. It so happens that the vacuum 
structure of the theory is that of a sequence of degenerate minima corresponding to integer 
values of the Chern-Simons number. Since these vacua are "pure gauge" , in the minima 
Chern-Simons number coincides with the topological winding number of the Higgs field; but 
whereas Higgs winding is always integer, Chern-Simons number can change continuously 
between vacua. 

Adjacent vacua are separated by potential barriers, and in equilibrium, the rate of 
transition from one to the other is a diffusive process governed by the Boltzmann factor of 
the saddle point configuration at the top of the effective potential, known as the Sphaleron 
[21, 22]. The Sphaleron energy is in turn dominated by the need for the Higgs winding 
number to wind in such a transition, since this can only happen through a local (in space) 
zero of the Higgs field length. The Sphaleron rate is by now very well known numerically 
[23, 24]. 



^This is in fact very similar to the situation in the Standard Model, except that there both the CP- 
breaking and the C-/P-breaking are associated with the fermions. 



Out of equilibrium, on the other hand, there is no concept of temperature and no Boltz- 
mann distribution. The dynamics of Chern-Simons number change need not be diffusive, 
and the change of Higgs winding number can happen spontaneously where the Higgs field 
happens to go through zero. This can happen through analogues of the Kibble mechanism 
[9, 25, 26], through overall oscillations of the Higgs field [8, 11, 27] or through statistical 
fluctuations as sphaleron-like half-knots [27]. In particular, Chern-Simons number and 
Higgs winding need not follow each other until the system begins to relax to equilibrium. 

In the 2HDM, we now have two Higgs winding numbers (one for each doublet), both 
of which must agree with the Chern-Simons number in the vacuum. Therefore, both must 
change in order to generate baryons, and so both Higgs fields must locally go through zero 
length during the evolution. This need not happen at the same time, nor at the same time 
as Chern-Simons number transitions, as long as they all end up at the same integer value. 
All in all this provides for some highly non-trivial dynamics. 

Spinodal transition: The Higgs field potential is such that at low temperature the elec- 
troweak symmetry is spontaneously broken through the appearance of a Higgs expectation 
value. The assumption is that at high temperatures this symmetry is restored and the 
expectation value vanishes as the effective potential changes. In the Standard Model we 
know that the transition is a cross-over [3], and that the rate of change of temperature in 
the Universe is governed by the Hubble rate. Since at the electroweak scale H ~ 10~^ eV, 
the system will remain very close to equilibrium throughout. 

Extensions of the scalar sector of the theory allows for a strongly first order phase 
transition, in which case bubble nucleation takes place, a very much out of equilibrium 
process. Baryogenesis can then occur at the bubble wall as it sweeps through the plasma 



We will be interested in another possibility, namely where the tree level potential is 
such that electroweak symmetry breaking is triggered from a super-cooled state. Consider 
as a toy model the potential 

V{ct), a) = {g^a^ - fi^) cj^U + A(<AV)' + ^o + V{a), (1.2) 

{V(a) is some appropriate potential for a only, and Vq fixes the minimum to have V = 0.) 
As long as {a) > fi/g, there is no symmetry breaking, and the Universe will cool down 
under expansion, in principle to arbitrarily low temperature. If the dynamics of a is then 
such that it goes from a > fx/g to cr ^ fj,/g very fast (on the time scale of /i), the <j) field will 
experience a fast mass quench and IR modes will become unstable. This can be seen from 
the linearized Klein-Gordon equations with negative mass-squared, in momentum space 



df + k^-i^^ 



^ 0k oc OkcV^'-'^'* + aj^e-Vfc'-M'*, (1.3) 



so that for |k| < /.i, occupation numbers grow exponentially 



nk oc e'V^^^^\ (1.4) 



until nonlinear interactions kick in and curb the growth. Large occupation numbers lead 
to approximate classical behaviour, which is why it is appropriate to evolve the system 
classically. A more detailed exposition of this point can be found in [29, 30]. Specific 
realizations of V{a) can be found in [31-33]. 

Classical Bosonic simulations: Very recently, an approach to simulate the whole 
SU(2)-Higgs-fermion system numerically has been developed [34-36], but the numerical 
implementation of the quantum fermions is very challenging in terms of computer re- 
sources. So for broad numerical investigations, it is vastly more efficient to consider an 
effective bosonic treatment, where the effect of fermions is captured through a series of 
higher order bosonic operators. 

As is always the case for effective bosonic simulations of electroweak baryogenesis, it is 
then assumed that when the gauge field Chern-Simons number changes in time, the baryon 
number of the fermions coupled to it follows along according to the anomaly equation (1.1). 
Hence we are interested in numerically computing {Ncs{t)), in our case using classical evo- 
lution equations for the fully coupled and non-perturbative SU(2)-2Higgs system, including 
CP-violation (from the tree-level scalar potential) and C- and P-violation (from effective 
higher order bosonic operators). 

The numerical procedure is then that the full quantum average of a certain operator 
is replaced by an ensemble average over field realizations 

l-^quant. ^ L^clas.ens. ; (-'^•"j 

and that each member of the ensemble is evolved independently using classical evolution 
equations. As mentioned, classical behaviour is applicable for fields with large occupation 
numbers as after a spinodal transition. 

The paper is organized as follows: In section 2 we will set up the 2HDM and the 
choice of scalar potential we use. In section 3 we describe the numerical setup and the 
initial conditions. In section 4 we present our numerical results, in 4.1 demonstrating why 
only CP-violation or only P- and C-violation does not generate an asymmetry and in 4.2 
performing simulations of a whole ensemble of configuration with both CP- and C- and P- 
breaking, computing the size of the asymmetry as a function of the strength of P-violation. 
We conclude in section 5. Some details of the choice of parameters for the Higgs potential 
can be found in appendix A. 

2 The 2HDM 

The 2HDM is defined through the continuum action 

S = - I d^xdt 



f 



^TrF^.^F'''^ + (L)^,,/.i)tZ)'Vi + (D^</)2)t 1)^02 + V{cPucP2) + Sc/p 



(2.1) 



where we use the metric rj = diag( — \- ++), 0i,2 are SU(2) doublets with hypercharge 
H-l and F^'^ is the field strength tensor of the gauge field. The covariant derivative is 



Dfj,<j)i = {dfj_ + iA^)(j)i and the potential is in all generality 

V {(1)1, (1)2) = -f4l4>\(t>l - f^22^2(t>2 - f^U <Pl4>2 - /^?2*4^1 

+ y ('AI</'i)' + y(4'/'2)' + A3(0l</.l)(402) + X,i4'PlM<P2) 

+y ('/'I<^2)' + ^{(P^lf + X6i4^i){4'p2) + A^(</.I</'i)(40i) 
+Xr{44)2){4<t>2) + X;{4<p2){4<Pi)- (2-2) 

The parameters Ai^2,3,4 and nfi 22 ^^^ ^^^^ ^^^'^ i^ general X^^j and /i^2 are complex. 

For the bosonic theory as written here, a generalized CP-symmetry can be defined 
which includes Higgs field basis transformations (in (/>i^2-space) [37, 38]. In this general 
sense, CP is then only violated if there is no such transformation that can make all the 
parameters real^. 

However, when the scalar fields are coupled to fermions in the usual way and the 
Yukawa couplings and CKM mixing matrix is fixed, there is no such freedom in general. 
In our case, we do not explicitly have fermions but imagine that they have been integrated 
out in the path integral, and that the leading effective bosonic operator resulting from this 
is [25, 26]3 

^c/P = T7r^^^(<^i'^2 - ^-^i) TrF^.F'^^ (2.3) 

iovr m^y 

where the Yukawa couplings and the mixing matrix is encoded in the real parameter (5c /p. 
This term breaks C and P maximally, but conserves the combination CP. Many more 
operators appear upon integrating out the fermions, but we will neglect all the C-even, 
P-even ones since they do not participate in generating an asymmetry, and take (2.3) 
as representative of the potentially many C- and P-violating terms at leading order. A 
complete calculation can be performed along the lines of [20], also providing a value for 
the coefficient S^/p. We again emphasize that we ignore the Standard Model CKM matrix 
as a source of CP violation. Had we included it, additional C-breaking, P-conserving (and 
hence CP-violating) terms would appear [20], but at higher order in gradients. 

From now on, we will set Xqj = for simplicity and we will fix the rest of the parameters 
in the following way (see also appendix A): 

• The five Higgs particle modes have masses m± = 400 GeV, rrii = 125 GeV, m,2 = 300 
GeV, 7713 = 350 GeV. The lightest mode represents the Standard Model-like Higgs 
[39, 40], and the remaining masses are chosen large enough to have avoided detection, 
while being "generic" in the sense of non-degenerate and not having some accidental 
symmetry. 

• The Higgs vevs are set by |i'i|/|f2| = tan/3 = 2 and vf + V2= 246^ GeV^. 



^We will not consider the possibility of spontaneous CP-symmetry breaking in this work. 
^Turok and Zadrozny [25, 26] consider a linearized, gauge fixed version of this term, where 
40i) ^ hihze, with 01 = (O./iie'"!)^, 02 = (0,/i2e^*2)^ and 9 ^ 61 -62. 



• The electroweak transition is driven by the spinodal instabihty of one or both of 
the Higgs fields. We will consider only the case when both of the eigenvalues of the 
matrix 

M'=(~^h~^lA (2.4) 

are negative. We indeed found that having both fields unstable (rather than only 
one) makes the system evolve faster and that the asymmetry becomes larger. 

• CP-conserving, parameter set 1: If all the parameters are chosen real, there is no 
CP-violation, and the vevs are real^. This leaves two free parameters. 

• Real vev, parameter set 2: We can have CP-violating also keeping both vevs real. In 
this case A5 is complex, but obeys the constraint Im(A5)wiW2 = Im(;Li^2)- We are left 
with three free parameters. 

• Complex vev, parameter set 3: Finally, we can make Higgs basis transformations and 
make A5 real, and the i?i)2-vev (for instance) complex. /i|2 is then also transformed. 
We again have three free parameters, and the potential is the same as for the real-vev 
case, but centered around the complex minimum. Then keeping (5c /p unchanged (the 
coupling to fermions), it is a different system to parameter set 2. 

We note that it really is C that is broken by complex parameters, whereas P is always 
conserved by the scalar potential. 

3 Numerical setup and initial conditions 

The continuum theory is discretized on a space-time lattice of volume L^ = a^n^ in a 
similar way to [11]. The lattice spacing is set through axV = 0.6, where v = 246 GeV, and 
the lattice size is Lv = UxttxV = 38.4, i.e. Ux = 64. In terms of the smallest and largest 
mass scales present, we have axini ~ 0.3, am± ~ 0.98 and Lrrii ~ 19.2, Lm± ~ 62.4 which 
we consider acceptable for the level of precision aimed at here^. 

The classical equations of motion are derived by straightforward variation of the action, 
and we consider the following observables: 



The Higgs field expectation values, normalized to their vacuum values 

_ 2 _ 2 

4>l = jY^^al(l)\(l)i{x), 4>l = jY^^ai(t>2<t^2{x), (3.1) 

'- X ^ X 



We also need a measure for the angle between the fields, which encodes the C(P) 
violation in time, as well as being the first factor of (2.3) 

lui{(t>l(t>i) = -^3^ Y. "' ^ ('^i'^^ - 4'^i) (^)- (3.2) 



*See footnote 2. 

^For classical simulations there is no need to take a strict continuum limit, as long as results vary little 
with lattice spacing and volume. 



In addition to the average, we will consider the distribution of the local quantity 
inside the sum. 



The Chern-Simons number 

iVcs(t)-A^cs(0) 



dt I d^x-^TrF^^F'"'. 



167r2 



fJ.U-' 



The Higgs winding numbers 



N" 



1 



247r2 



d xeijkTr 



UldiUaUldjUaU^adkUa 



a=l,2, 



where 



$. 



i>a,io-2(t)l), Ua 



$n 



:^n 



(3.3) 



(3.4) 



(3.5) 



The lattice winding numbers are integer up to discretization errors and both coincide 
with the Chern-Simons number in the vacua. 

• As a check on the numerical code, we keep track of Gauss constraint and energy 
conservation. 

We generate a set of random initial field configurations in the following way: Given 
the masses and couplings, we diagonalize the Higgs mass matrix (2.4), to find 



Mj 



Ml \ 
Mil 



(3.6) 



in terms of the eigenbasis 4>a,b- We restrict ourselves to cases where M^ and Af^ are 
both negative, and so both eigenmodes will have spinodally unstable (lattice) momentum 
modes, given by^. 



^lat - 2^ „2 



, 27r 
2 — 2 cos I —rii 



<\M\nl i = 1,2,3, ni = 0,...,n^-l. (3.7) 



These (and only these) momentum modes are then initialized in the "quantum vacuum", 
with a Gaussian distribution according to [29, 30] 



1 



{<t^A,B{k)^\Bik)) = l, 



1 



'^Jkh + \Ml 



{dtcl>AAk)dtAA^)) = \^kl^ + W. 



2 I 

a,b\ 



"lat 



'A,B\ 



This mimics a fast mass quench, where the potential instantaneously flips as 

t<0, n<^4,<^B) = ((^^.,(/.Ul '^^' ,,°., 1 f t^ V (3-8) 



WP 



B\ 



t >0, 



V{ct>AAB) = {ct>U^s)\^^^ ^p^ 



+ quartic. 



'm 



(3.9) 



"if only one of the eigenvalues M\ g is negative, only the unstable modes of the corresponding field 
should be initialized. 



The initialized fields are then rotated back to the original i;^i^2 basis. 

This is completely analogous to the approach in [11], generalized to more than one 
Higgs field. Similarly, the gauge field is initialized Ai = 0, and the simulations are done 
in temporal gauge Aq = 0. The gauge momenta dfAi are initialized by solving Gauss 
constraint in the background of the random Higgs field. 

Since we are interested in C-, P- and CP-violation, and in order to minimize the 
statistical noise, we generate an ensemble of initial configuration that is C-even and P- 
even. This is of course automatically satisfied for an infinitely large ensemble, since a 
configuration and its C- and/or P- conjugate have the same probability to be generated. 
But we make it explicit by, for all random configurations (/)init(x), ^init(x) (in terms of 
link variables [/j(x)), also including their C-, P- and CP- conjugate configurations in the 
ensemble, according to 

$i(i,f) A$,(t,-x), Un{t,x)^Ul{t,-x-h), (3.10) 

and 

^i{t,x)^^*{t,x), Un^U*{t,x). (3.11) 

In the continuum limit, this reduces to 

$i (t, f ) A $i (t, -f ) , A{t,x)^-A{t,-x), (3.12) 



and 



$i {t, x) A «^* (t, x) , A^ {t, x) A -A* {t, x) . (3.13) 



CP-conjugation is of course the combination of the two. We then construct averages for a 
given observable O through 

= ^^0{{^,}), (3.14) 

-'Vinit 

I 

where 

O({0a) = \ [0{ct>i) + Oic^f) + O(0f ) + O(0f ^)] . (3.15) 

This is identical to a normal average over all 4x A'^init configurations. But for the statistical 
error we use 



o^im) - io{{cp}) 

4 = ^7^ '-■ (3.16) 

We can convert the non-zero average Chern-Simons number to a baryon asymmetry 
using the conversion prescription (we use the Standard Model content of fermionic and 
bosonic degrees of freedom) , 

^=(2^wfe- P^' = nm-Vin,..), (3.17) 
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Figure 1: The Chern-Simons numbers, winding numbers and Higgs expectation values 
in time for a single configuration. Left: CP-conserving potential. Right: CP-breaking. 



i.e. distributing the initial potential energy on a thermal ensemble of all the relativistic 
degrees of freedom g* . For the Standard Model degrees of freedom with masses less than 



m^ 



we have g* = 86.25. This amounts to 






X 1.2 X 10" 



V{^,{))-V{VI,V2) 



-3/4 



The reheating temperature is, for the three choices of potentials respectively, 

T 



GeV 



(99, 104, 104) . 



(3.18) 



(3.19) 



This is much higher than in the Standard Model, simply because the potential is deeper 
(we may even consider treating the W^ and Z as relativistic, changing the temperature 
by a factor (95.25/86.25)^'^ = 1.025, i.e. 2 — 3 GeV lower temperature). Consequently, 
the thermal electroweak transition is also at higher temperature than in the Standard 
Model, and the sphaleron rate is modified accordingly. A detailed study of these features 
of the 2HDM model would certainly be of great interest both for hot and cold electroweak 
baryogenesis. 

There is an interesting twist, in that for the simulations performed here, the total 
number of degrees of freedom participating in the dynamics is g* = 14 (or g* = 10, if we 
neglect the heavy Higgs modes), in which case the effective temperature goes up by a factor 
(86.25/14)^'^ ~ 1.6, which brings it above the lowest Higgs mass of 125 GeV. This means 
that the temperature is large enough to partially "restore" symmetry in one direction in 
field space. This is an (albeit physical) artefact of not including all the existing degrees of 
freedom in the dynamics. We will comment further on this below. 



4 Results: 



4.1 



Single trajectories and symmetries 

We will start out by considering a single initial field realization, and in Fig. 1 (left), we 
show the Chern-Simons number Eq. (3.3), both Higgs winding numbers Eq. (3.4) and both 



Higgs expectation values Eq. (3.1). We use the CP-conserving potential (parameter set 1). 
We observe that both Higgs fields perform a spinodal transition in a time t ~ W/v, ending 
up close to their respective vevs (here normalized to 1). Chern-Simons number oscillates 
and drifts in to A'^cs — —2, and around t = 20/w, both Higgs winding numbers settle close 
to Chern-Simons number. We also see that N^ is still quite noisy at later times, while 
N^ stays flat. This signals that the mass in the "1" direction is smaller than in the "2" 
direction, and that it is therefore easier for the field to make excursions near (^i(a;) ~ 0, 
where the lattice winding numbers are ill defined and easily "flips" . The winding numbers 
a slightly off integer values because of discretization errors; the Chern-Simon number is 
not exactly integer because the fields have finite energy density/temperature. 

In the right-hand figure we show the same observables, but for another single config- 
uration evolving in a CP-breaking potential (parameter set 2). There is now a substantial 
difference in that the "1" mode is much lighter than the "2" mode. The "2" mode is 
strongly forced into its minimum, whereas the "1" field has large amplitude oscillations. 
The corresponding winding numbers are also different, in that A^^ is smooth and closely 
integer, whereas N^ is noisy and has spikes and nearly jumps. This is because the field 
oscillation brings about zeros of (j)i{x). The spikes are also strongly correlated with the 
overall oscillations of the field. 

In Fig. 2 we show the Chern-Simons number and winding number "2" for four conjugate 
initial configurations (base, its C-, P- and CP-conjugate), the P- and CP-conjugates with 
an overall sign-flip. In the top left plot, the potential conserves CP and the C-/P- breaking 
term is turned off (^c/P — 0- The curves are identical, and we conclude that the average 
Chern-Simons number is identically zero, since the P- and CP-conjugate exactly cancel the 
other two. In other words, iVcs({</'}) = N^iW) = 0. 

In the top right plot we show the Chern-Simons and winding numbers in a similar run, 
but now in a CP-breaking, real-vev potential (parameter set 2). The four configurations 
are no longer identical up to the overall sign, but pair up in P-conjugates. The base 
configuration cancels with its P-conjugate; the C-conjugate pairs up with the CP-conjugate. 
Again the averages Ncs{{(l)}) and N^{{(j)}) are zero. 

Then in the bottom left plot we have included the P- and C- violating term Eq. (2.3), 
but the potential is again CP-conserving. Now the pairing is the other way around: The 
base configuration is cancelled by the CP-conjugate, whereas the C- and P-conjugates 
average out. Ncsdcp}) = N^{{(f)}) = 0. Finally, including both the CP-breaking potential 
and the C-/P-breaking term, we find as shown in the bottom right plot that the four 
trajectories of Chern-Simons number are all different. Ncs{{(l>}) 7^ and N^{{(f)}) ^ 0. 
Hence a net baryon asymmetry is created. 

4.2 Ensemble averages 

In order to compute the expectation value of the baryon asymmetry, we now perform sim- 
ulations of an ensemble of 4 x 25 configuration (4 for the conjugates) and compute average 
Chern-Simons number and Higgs winding number, as shown in Fig. 3 for the complex- 
vev (left) and the real-vev (right) CP-breaking potentials. The error bars/band on N^ 
correspond to standard deviations as defined in Eq. (3.16). 
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Figure 2: The Chern-Simons number and winding number "2" for a single configuration 
and its C-, P- and CP- conjugate in tlie case of no C-, P- or CP- violation (top left) 
and when including CP-violation in the scalar potential (top right). Also the case of C-, 
and P-violation but no CP-violation (bottom left) and when including CP-violation as 
well (bottom right). The P- and CP- conjugates have been given an overall sign flip for 
illustration. 



For the complex-vev case, we see a large initial creation of Chern-Simons number, 
driven by the C-/P- violating term as the Higgs fields go through the transition. This 
initial growth is not matched initially by the winding numbers, but around time t = 20/t;, 
these also acquire a non-zero average value. Eventually by t = 30/t;, all three observables 
have settled to a common asymmetry, and nothing more happens. The only remaining 
effect is that Chern-Simon number oscillates in time with the Higgs fields, in effect the 
coefficient of the C-/P-breaking term. 

For the real-vev case, there is no initial bouncing of the Chern-Simons number, and 
together with the winding numbers there is a slow drift to a finite value, which within error 
bars matches the complex-vev result. This may seem surprising, but means that whether 
or not the late-time Higgs expectation value is real or complex is not very important for the 
asymmetry. Baryogenesis takes place in the initial very chaotic stage, where CP-violation 
is active through complex couplings irrespective of the reality of the vevs. 

P violation results from the C-/P-breaking term when the coefficient Im((/)2(^i) = 
I (^1 1 1 (^2 1 sin is non-zero, effectively biasing Chern-Simon number. In Fig. 4 (left), we 
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Figure 3: Chern-Sinions number and Higgs winding numbers when averaged over an 
ensemble of 4 x 25 configurations. The potential is CP-breaking with complex vevs (left) 
and real vevs (right). 5q/p = —105. We have included the la statistical error band for N^ 
for illustration. The errors on the other observables are similar. 
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Figure 4: Left: The evolution of the normalized Higgs expectation values and the av- 
erage complex angle 6 for a CP-breaking, and complex-vev potential, and its evolution in 
time. Right: The distribution of the Higgs field imaginary part |(^i||(^2| sin 6 at the initial, 
intermediate and final times (t =0, 2, 5, 120). 



show the evolution of the (normalized) average Higgs fields and the relative angle 9 for 
the complex-vev potential. All settle fairly early on, t ~ W/v, resulting in a non-zero P- 
violating coefficient. In Fig. 4 (right) we show the distribution of Im((/>2</>i)(x) at the initial, 
final and two intermediate time-slices. We see how the initial condition is C-symmetric and 
strongly peaked, and then as the system evolves, the distribution fiattens out and moves 
to a non-zero average value. 

For the real-vev case, CP-violation shows up as a nonzero average of the Higgs field 
angle Eq. (3.2) at intermediate times. In Fig. 5 (left) we show this average, and in Fig. 5 
(right) the distribution at different times. The distribution is again symmetric and peaked 
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Figure 5: Left: The average complex angle for a CP-breaking, real-vev potential, and 
its evolution in time. Right: The distribution of the Higgs field imaginary part |0i||02| sm.9 
at the initial, intermediate and final times {t =0, 2, 5, 120). 



initially, but now during evolution the average moves away from zero, and then returns. 
The distribution again flattens. We also see that the late times average value is no longer 
exactly zero, showing that the minimum of the effective potential at finite temperature 
moves to complex expectation values, even though at tree level the vev is tuned to be real. 
In Fig. 5 (left) we also shos the average fields, and again that the potential is rather shallow 
in the (/)i-direction, so that at finite temperature the field oscillates with large amplitude. 
The final baryon asymmetry depends on the dynamics of the system, and hence the 
precise form of the potential. For the present work we will not sweep through this rather 
large parameter space, but instead fix the potential, and calculate the dependence of the 
asymmetry on the parameter (5c /p. The net Chern-Simons and winding numbers are shown 
for the complex- and real-vev cases in Fig. 6. We first remark that because at late times the 
Chern-Simons number is still oscillating, the cleaner observables are the winding numbers. 
These have a smooth, monotonic but non- linear behaviour as a function of ^^/p, which 
goes very neatly through the origin as it must by construction. 

— 1 2 

Based on the result for N^ at ^c/P = ~21, we conclude that the baryon asymmetry 
is given by 






for the complex-vev potential and 

riB _ 



-(5c/p X (2.8 ± 1.2) X 10^^ 



-5c/p X (1.6 ±1.2) X 10" 



(4.1) 



(4.2) 



for the real-vev potential. These two results are consistent with each other at the la level, 
and this agreement applies to the full (5c/p-dependence in Fig. 6. Again, it seems to have 
little impact whether in an otherwise identical potential, the expectation value is rotated 
to a complex value or not. 
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Figure 6: The dependence of the baryon asymmetry on the strength of P-violation, for 
the CP-breaking, complex vev (left) and real vev (right) potentials. 



5 Conclusion and outlook 

Using fully non-perturbative and out-of-equilibrium classical lattice simulations, we have 
computed the baryon asymmetry generated in a cold spinodal electroweak transition in 
the Two-Higgs-Doublet model. CP-violation appears as C-violation in the scalar potential, 
but only in combination with explicit P (and in our case, also C) breaking can a net 
Chern-Simon number density be created. In the full theory, this is present through the 
usual left-handed projector in the gauge-fermion interaction. In this bosonic model, this 
is replaced by the leading C-/P-breaking bosonic operator expected upon integrating out 
the fermions. 

Our results show that an asymmetry is indeed created; it is large enough to be seen 
clearly in simulations, and is not suppressed by being a combination of two (small) effects, 
the CP-violation and the C-/P-violation. The observed cosmological baryon asymmetry of 
~ 6 X lO^'^'^ is reproduced for 



5c/p = -(2 to 3) X 10-^ 



(5.1) 



The asymmetry is largely created in the initial spinodal roll-off, and the late time dynamics 
is mostly irrelevant. In particular, whether or not the Higgs vev is real or complex had 
little impact (if the potential is otherwise the same) . 

In the potentials used here, the temperature after the transition is high enough that 
there is "partial symmetry restoration" in the lightest Higgs mode direction. This is partly 
because all the fermion degrees of freedom are not included dynamically and g* is much 
lower in the simulations. This restoration makes it possible for one of the Higgs fields 
to unwind, but since both winding numbers and Chern-Simons number are dynamically 
chained to each other at late times, this does not happen. At early times, Chern-Simons 
number and the winding numbers move independently. 

The "missing link" is how to trigger the required spinodal transition and how the 
Universe came to be super-cooled way below the electroweak scale. This can be achieved 
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Figure 7: The potential for parameter set 1, CP-conserving, real vevs. Plotted in ui,2- 
space at ^2 — ^1 = 0. 



through a further enlarged scalar sector, for instance by a gauge singlet, which may [32] 
or may not [31] be taken to be the inflaton. Further examples of potentials, which may 
incorporate a first order transition as well, have been proposed in [33]. 

The obvious next step is to compute the value of ^c/p, or in more generality, the 
functional form and coefficients of all C- and P- violating terms arising from integrating 
out the fermions. This can for instance be done using a gradient expansion as in [14- 
17, 20]. The precise asymmetry generated will depend on the exact scalar potential, but 
we expect the Sqjp dependence found here to persist for a given potential. A more detailed 
study of the space of scalar potentials, including cases with only one unstable initial field, 
would also be a natural extension of the present work. Finally, the inclusion of the U(l) 
gauge field, a dynamical scalar to trigger the spinodal transition and ultimately dynamical 
fermions would all be interesting extensions and improvements. 
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A Choice of Higgs potential parameters 

We consider three different Higgs potentials, denoted as parameter sets 1, 2 and 3 respec- 
tively. They all obey our main criteria for the parameters Ai^2,3,4,5 and the mass parameters 
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Figure 8: The potential for parameter sets 2 and 3, CP-breaking, real/complex vevs. 
Plotted in vi 2-space, at ^2 — ^i = and 1.39, respectively. 



Mil 22 12 that only /ifg ^^'^ -^5 ^^^ be complex. We impose the following constraints (See 
also [37, 41, 42]): 

• The global minimum has (j)i = (0,fi/\/2e*^^)"^, (f)2 = (0,f2/\/2e*^^)"^, with Ui/f2 = 
tan/3 = 2 and ^^vf + v'^ = 246 GeV, so that vi = 220 GeV, ^2 = HO GeV. 

• The mass eigenvalues around the global minimum are (m±,mi,m2,m3) = (400, 125, 300, 350) 
GeV. 

• In the CP-conserving case, both A5 and /_f^2 ^^^ real, and 9i — 62 = 0. 

• In the real-vev case we choose Oi — 62 = 0, in which case we have the constraint 
ImA5fit;2 = Ini/if2- 

• In the complex-vev case 9i — 62 is free and we choose A5 to be real. 

This leaves two, three and three parameters, respectively to be fixed, and we use the sets 
(masses squared in GeV^) : 

• 1) CP conserving: (see Fig. 7) 

Ai = 2.5575, A2 = 1.5424, A3 = 6.1176, A4 = -3.0570, A5 = -1.8177, 

nl^ = 71125, HI2 = 84713, ^^2 = 10000 (A.l) 

The eigenvalues of the initial condition mass matrix are (M^, M^) = (300.0^, 256.6^) 
GeV^ for which the lattice size of (Lv)^ = (38.4)^ is easily large enough for sufficiently 
many spinodal modes to be present. 
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2) Real-vev: (see Fig. 8) 

Ai = 0.86175, A2 = 2.36749, A3 = 5.7886, A4 = -3.5845, A5 = -1.2902 + i0.48231, 
^2^ = 34673, /i^2 = 120680, fij^ = 10000 + ill675 (A.2) 

The eigenvalues of the initial condition mass matrix are (M^, M^) = (351.2^, 178.9^) 
GeV^ for which the lattice volume is again large enough. 

3) Complex-vev: (see Fig. 8) 

Ai = 0.86175, A2 = 2.36749, A3 = 5.7886, A4 = -3.5845, A5 = 1.3774, 

^2^ = 34673, 1^4^ = 120680, fij^ = 13268 - i7763. (A.3) 

The global minimum has 62 — 9i = 1.39. The third potential is simply a rotation 
of (j)2 by an angle 82 — Oi to make A5 real. So the potential and its minimum are 
equivalent, if the P-breaking coefficient were also rotated by the same complex phase. 
If Sq/p is kept constant, it is a different physical system. 



f V{0,0)-V{v„V2y , ^^^^^^^ Q_g^Q^ Q g^Q^ ^^ ^^ 



We have for the three sets, respectively 

3/4 

V yu, Kjj — V \^L/i, U2) y 

which enters in the calculation of the final asymmetry Eq.(3.18). 
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